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Abstract: 

We consider quantum interferences of classically allowed or forbidden electronic tra- 
jectories in disordered dielectrics. Without assuming a directed path approximation, we 
represent a strongly disordered elastic scatterer by its transmission matrix t. We recall 
how the eigenvalue distribution of t.t^ can be obtained from a certain ansatz leading to 
a Coulomb gas analogy at a temperature (3^^ which depends on the system symmetries. 
We recall the consequences of this random matrix theory for quasi-ld insulators and we 
extend our study to microscopic three dimensional models in the presence of transverse 
localization. For cubes of size L, we find two regimes for the spectra of t.t^ as a function 
of the localization length C,. For L/^ ^ 1 — 5, the eigenvalue spacing distribution remains 
close to the Wigner surmise (eigenvalue repulsion). The usual orthogonal-unitary cross- 
over is observed for large magnetic field change AB ^ $o/^^ where $o denotes the flux 
quantum. This field reduces the conductance fluctuations and the average log-conductance 
(increase of ^) and induces on a given sample large magneto- conductance fluctuations of 
typical magnitude similar to the sample to sample fluctuations (ergodic behaviour). When 
^ is of the order of the lattice spacing (L/S, ^ 5 — the eigenvalue repulsion is weaker 
and the removal of a time reversal symmetry has a more negligible role. In those two 
regimes, the sample to sample fluctuations of the log- conductance are close to a normal 
one-parameter distribution, and depends linearly on the disorder parameter. In these 
microscopic models, the Fermi energy dependence of the conductance is very similar to 
the one recently observed in small GaAs:Si wires. 
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Resume: 



Nous considerons les interferences quantiques de trajectoires electroniques classique- 
ment permises ou interdites dans des dielectriques desordonnes. Sans faire une approxi- 
mation de chemins diriges, nous representons un diff'useur elastique tres desordonne par 
sa matrice de transmission t. Nous rappelons comment la distribution des valeurs propres 
de t.t^ peut etre obtenue d partir d'un certain ansatz conduisant a une analogic avec un 
gaz de Coulomb a une temperature qui depend des symetries du systeme. Nous rap- 
pelons les consequences de cette theorie de matrices aleatoires pour des isolants quasi-ld 
et nous etendons notre etude d des modeles micro scopiques tridimensionnels ou la locali- 
sation transverse est presente. Pour des cubes de taille L, suivant la valeur de la longueur 
de localisation ^, nous trouvons deux regimes pour les spectres de t.t'^ . Quand L/^ ^ 1 — 5, 
la distibution des espacements entre valeurs propres reste proche de celle de Wigner (re- 
pulsion des valeurs propres). Le cross-over habituel entre les cas orthogonal et unitaire 
a lieu pour un champ magnetique applique eleve AB ^ ^o/^^ ou $0 est le quantum, de 
flux. Ce meme champ reduit les fluctuations de conductance g ainsi que la moyenne de 
log{g) (augmentation de ^) et induit sur un echantillon donne de grandes fluctuations 
de la magnetoconductance d'une ampleur comparable a celles d' echantillon a echantillon 
(comportement ergodique). Quand ^ est de I'ordre de la maille du reseau, (L/^ » 5 — Gj, 
la repulsion des valeurs propres est plus faible et I'effet du champ magnetique cesse d'etre 
important. Dans ces deux regimes, les fluctuations d'echantillon a echantillon de log{g) 
sont bien decrites par une gaussienne a un parametre et varie lineairement en fonc- 
tion du parametre de desordre. Dans ces ces modeles microscopiques, la dependance en 
fonction de Venergie de Fermi de la conductance est tres comparable a celle recemment 
observee sur de petits fils de GaAs:Si. 
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1. Introduction 

When one has to tackle the problem of electronic quantum transport, one of the usual 
microscopic starting points is the familiar Anderson tight-binding Hamiltonian: 



where the are the random site energies and the hij represent the nearest neighbor 
couplings or transfer terms. For instance, the can be random with a rectangular distri- 
bution of width W (Anderson model) or can take at random two values or ej, (binary 
alloy model). The h^j can take for simplicity a uniform value h ii i,j are nearest neigh- 
bors, otherwise. Perturbations could be done when one of the two terms of H (random 
substrate electrostatic potential or kinetic energy) is neghgible compared to the other. 
The scaling theory of localization [1] assumes the existence of a universal /9(5')-function 
of conductance g relating those two limits through the mobility edge in three dimensions. 

For disordered metals whose quantum conductance g, measured in unit of e^/h, is 
larger than one, the conventional perturbation approach starts from plane-wave like 
states which are only weakly perturbed by the scattering caused by the random substrate 
(W/h ^ 1). In this limit, electronic transport essentially results from quantum interfer- 
ences of classically allowed diffusive paths and quantum theory is close to its semi-classical 
limit. 

(i) We know for the three different possible symmetries (orthogonal, unitary and sym- 
plectic cases): the four loop order expansions [2] in 2-|-e dimensions of the scaling functions 
P{g). One gets after integration the weak localization corrections to Boltzmann conduc- 
tance, i.e. how < g > depends on the sample size L, the brackets denoting the ensemble 
average value. 

(ii) The sample to sample conductance fluctuations around the ensemble average essen- 
tially satisfy gaussian distributions whose variance depends only on the basic symmetries 
and slightly on the dimensionality of the system [3] (Universal Conductance Fluctuations). 

(iii) The conductance fluctuations induced on a given sample by the variation of an ap- 
plied magnetic field or of the Fermi enegy are of similar magnitude [4] (crgodic behaviour). 
An applied magnetic field removes time reversal symmetry, suppressing the weak (anti) 
localization correction for < > in the absence (in the presence) of spin-orbit scattering, 
halving the variance of the universal conductance fiuctuations. The cross-over field is in 
general small in metals, corresponding to the application of a fiux quantum $0 through 
the temperature dependent area where electrons keep their quantum coherence. 

(iv) The conductance fluctuation yielded by a local reorganization of the substrate 
potential (single impurity move) have been understood mostly from known results on 
random walks, yielding a possible explanation for universal amplitude of observed 1//- 
noises [5]. 

For Anderson insulators, W/ h^ 1 and it is better to perturb the point-like-localized 
eigenstates of the random term in (1) by the kinetic transfer term. The perturbation 
removes the possible degeneracies of the e, (which is huge in the binary alloy model) and 
the eigenstates of (1) are still localized over domains of typical size C, which can extend 
over many impurities in the vicinity of the mobility edge. The Green function G{E) 
between two sites i and j can be expanded in powers of the hopping integral h 
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Here we have also introduced a magnetic vector potential A for generality. The paths T 
going from i to j are in this case a priori highly non classical, and electronic transport 
mainly results from classically forbiden trajectories which interfere due to the random sign 
of the denominator and to the vector potential. When all the factors h/ {E — e^) <^ 1, the 
shortest directed paths Vd need only to be considered [6] and electron trajectories with 
self intersections have no significant contribution in this limit: 

lA 

<.iG(E)ij>^/.^EnT^- (3) 

^r, ^ - 

The combinatorics of the shortest directed paths of length L can be simply obtained 
if some cases (e.g if the sites i and j stand at the opposite corners of a square lattice). 
More generally, these directed paths relating two points separated by L have transverse 
fluctuations extending over {L^Y- Assuming standard diffusion, one gets ^ = 1/2 but a 
more detailed analysis shows that dominant paths are "superdiffusive" , yielding different 
values [7] for C,. One finds that: 

(i) The locahzation length ^ characterizing the decay of the propagator: 

< log\ < i\G{E)\3 > \^ >^ L{p - i^') (4) 

depends on the disorder parameter W through a local contribution = \!'Og{W / h)]~^ 
which comes from the exponentially small hopping term [h/WY" in (3) (binary alloy 
model with ej =1 W aX, E = Q). po is a global contribution embodying the M^-independent 
quantum interferences of the directed paths in this particular model [8]. ^ increases as 
a function of an applied field B as B^/'^ with a universal factor, giving rise to a posi- 
tive magnctoconductance [9]. This sign of the magnetoresistancc is not reversed by the 
presence of a strong spin-orbit scattering [10], as it is in disordered metals. 

(ii) The sample to sample conductance fluctuations are lognormally distributed, but 
with anomalous size dependence of the fluctuations. var{log{g)) = C.L'^'^, where the 
prefactor C does not depend on the mean < log{g) >oc L, a; = 1/3 in two dimensions and 
a; ~ 1/5 in three dimensions [7-8]. The distribution of log{g) then contains more than a 
single parameter. 

(iii) The application of a flux quantum $o through the "cigare-shape domain" , which 
contains the forward scatterered trajectories, yields a fluctuation ^(^f) g much smaller 
than the typical sample to sample fluctuations [11] (non ergodicity). No other character- 
istic magnetic field scale appears if self-intersecting trajectories are neglected. 

(iv) The effect of the motion of a single scatterer has not been addressed in directed 
path theories. Broad distributions of the log-conductance with a quasi-universal variance 
have been numerically observed without a directed path approximation [12]. 

The restriction of quantum interferences to this class of directed trajectories is usually 
made in ad hoc models (binary alloy distribution at an energy far from the e,,) where 
one assumes that the path contribution is always an exponentially decreasing function of 
its length. This does not take into account situations where small denominators could 
compensate the exponentially small value of h^. Such quantum resonances necessarily 
occur in a device where carrier concentration can be varied by capacitive technics. This 
approximation does not describe also more general cases where transport involves com- 
plicated interferences between classically forbidden and allowed trajectories (vicinity of 
the mobility edge). In an Anderson model with a broad distribution of h/ci, there is a 
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finite probability to have clusters where the electron can diffuse easily, be backscattered 
and selfintersects its own trajectory before having to tunnel through large energy barri- 
ers. Such metallic clusters can yield non trivial magnetic field or spin-orbit scattering 
dependence in the vicinity of the mobility edge. 

An alternative non perturbative theory for the analysis of quantum transmission in 
insulators is provided by a maximum entropy ansatz [13] using the transfer matrix M. 
This approach quantitatively describes quantum transmission through quasi-ld wires over 
length scale smaller or larger than their localization lengths, leading us to predict novel 
universal symmetry breaking effects [14]. We are interested to know its degree of validity 
outside quasi-one dimension in the presence of a strong transverse localization, and to 
compare its predictions from perturbation theory restricted to directed paths. We do 
not restrict ourself to the limit of an elementary localization of the conduction electrons 
around randomly located donors or acceptors, i.e. to a purely geometrical disorder, but 
we consider also systems closer to a mobility edge where localization results from complex 
interferences between classically allowed or forbiden trajectories. In this later case, it is 
likely that the relevant quantum interferences do not only involve simple forward scatter- 
ered directed paths, and that the effect of a strong spin-orbit coupling or of a large applied 
magnetic field have to be revisited [15, 16]. Being outside a simple perturbative limit, 
we rely on numerical simulations which include all possible path contributions without 
any a priori simplification. We use as microscopic models either a random quantum wire 
network or an Anderson model with rectangular distribution of ej. 

This work belongs to a series of three papers devoted to quantum interferences in dis- 
ordered insulators. One of them is theoretical [17] and analyzes how the purely quantum 
interferences which we consider are truncated by activated process at low temperature. 
The second [18] is experimental and gives the energy and magnetic field dependence of 
the conductance. One could then compare the results of our simulations to their results 
and note a striking (at least qualitative) similarity. The difficulty for a more quantitative 
comparison is due to the activated processes. 

We recall in sections 2 and 3 the maximum entropy description of quantum trans- 
mission and its implications for quasi-one dimensional insulators (typical conductance 
and sample to sample fiuctuations, symmetry breaking effects). We study in section 4 
the modification of the physical constraint used in this maximum entropy description for 
three dimensional insulators. The persistence of the short range repulsion in presence of 
transverse localization is shown in section 5 through a study of the eigenvalue nearest 
spacing distribution of t.t"^, though the long range interactions are weaker than usual 
random matrix behaviour (study of the spectral rigidity). In section 6, the sample to 
sample fiuctuations of ln{g) arc studied in three dimensions for large disorder. We show 
that varies as a function of the disorder parameter linearly from the mobility edge to 
the strongly localized regime. The effects of a variation of the applied magnetic field and 
of the Fermi energy are considered in section 7 and 8. We compare our results for the 
eigenvalues of t.t''" with those which are understood for the system Hamiltonian (section 
9), and with the experimental results of Ref. 18 (section 10). 

2. Mctximum entropy description of quantum scatterers 

We consider a A'"-channel quantum scatterer that we represent by its transfer matrix 
M. This matrix gives the electronic fiux amplitudes at the right hand side of the scatterer 



5 



in terms of the flux amplitudes standing at its left hand side. The statistical ensemble for 
a matrix X, related to M by the relation 

X = {M\M +{M.M^)-^ -1)/A, (5) 

is chosen [13] as the most random matrix ensemble given the density o-(A) of the eigenval- 
ues of X. This statistical ensemble of maximum information entropy for X implies that 
the X-eigenvectors are statistically independent of the eigenvalues and are distributed 
totally at random. We denote by t the transmission matrix at the Fermi energy. Us- 
ing a two-probe Landauer formula we express the conductance g — 2.ir(t.t^) as a linear 
statistics of the eigenvalues Ta of t.f^ or the eigenvalues {Aa} of X via the relations: 



N N 1 

g^'2.Y.Ta = 2.Y.-—-, (6) 

1 „_ i i -r An 



a=l a=l 



the factor two coming from a possible spin or Kramer degeneracy. Sometimes, we refer to 
the eigenvalues {Aa} of X as the radial parameters of M, since they can be also introduced 
through a natural "polar" decomposition of the transfer matrix. 

The joint probability distribution P({Aa}) can be easily calculated from this maximum 
entropy ensemble. One gets an expression formally identical to the Boltzmann weight of 
a set of classical charges free to move on a half line (the positive part of the real axis) at 
a well-defined "temperature" equal to 

P({A})«exp-/3.i/({AJ), (7) 

and governed by an hamiltonian: 

N N 

H{{\a}) = - E In |Aa - M + E V{\a). (8) 

a<h a=l 

The confining potential is given in terms of the input density (t(A) by 

roo 

V{\) = / cT(A')ln|A- A'ldA'. (9) 

Jo 

In usual random matrix theories, this is the classical Coulomb gas analogy which 
introduces a temperature, a one-body potential and a two-body interaction. 

(i) The temperature is given by the fundamental symmetries which leave the system 
invariant (time reversal symmetry, spin rotation symmetry), f3 = 1,2,4 respectively for 
the orthogonal, unitary and symplectic case. In the absence of spin-orbit scattering, 
(3=1 without applied magnetic field and becomes equal to 2 if the applied magnetic field 
is large enough (removal of time reversal symmetry). 

(ii) The confining potential V{\) can be adjusted in order to attract most of the Aa 
near the origin, yielding a large value for g (metal), or on the contrary V{\) can be very 
weakly attracting, leaving the "charges" spread out along the positive real axis, yielding 
g <^ 1 (insulator). In other words, the spectrum of X depends on the different physical 
parameters (Fermi momentum kj, disorder measured by system size L...) through a 
simple one-body potential V{\) in this maximum entropy ansatz. 

(iii) As the value of /3, the crucial logarithmic interaction between the A^ is obtained 
from the invariant measure of the space of matrices X which are invariant under certain 
symmetries that each member of the statistical ensemble must satisfy. 
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The consequences of the ansatz (8) on the spacing distribution of successive radial pa- 
rameters and on their two point correlation functions have been carefully investigated [19- 
20] for disordered conductors where the Aa are in part accumulated close to the origin. 
Detailed studies confirm both the presence of the infinite range logaritmic interactions 
and the change of temperature associated with symmetry breaking effects when the sam- 
ple shape is nearly quasi-one dimensional. For higher dimensions {d = 2 and 3), the 
transverse diffusion slightly reduces the rigidity of the A-spectra for high transmission 
or reflection (edges of the spectra), albeit the random matrix ansatz remains a good 
approximation. 

In the quasi-l(i limit, the ansatz (3) is consistent with a diffusion equation [13] which 
implies the quantitatively correct behaviours for the average and the variance of g in the 
metallic limit, if the input density cr(A) satisfles a certain integro-differential equation. 
More important for the questions investigated here, this ansatz remains valid when the 
sample length is larger than the quasi-ld localization length, establishing the validity of 
this "Coulomb gas model" in the presence of longitudinal localization. In other words, 
the Coulomb gas model remains essentially stable with respect to matrix multiplication 
for a flxed number N of channels. 



3. Longitudinal localization without transverse localization (quasi-ld limit) 

We summarize in this section previous results obtained from the global ansatz [13-21] 
and based on an analytic expression for (t(A) yielded from the diffusion equation [22] valid 
in the quasi-lo? limit. The analysis is easier in terms of the variables {va} related to the 
{Aa} through 

A. = (10) 

In the large A^-hmit, one gets for the i/-density 

21 

a(u) = 0<u < 2L/1 (11) 

and zero elsewhere which gives for the confining potential: 

U(iy) = / In 1 cosh(i/) - cosh(iy')\diy' ^ --^v^. (12) 
2L Jo 4L 

We write again the distribution P({z/(j}) as a Boltzmann weight at a temperature (3"^ of 
a classical system, getting for the hamiltonian 

TV N 

H{M) = - E M^a, l^b) + E ^l(^c). (13) 
a<b c=l 

The two body interaction is now: 

V2{va, J^b) = In ] cosh(i/a) - cosh(i/b)], (14) 
and the one body potential contains a /3-dependent term yielded by the change of variable: 

V,{iy,) = ln(sinh(z/e)) + U{u,). (15) 
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However, the quasi-lci localized limit is characterised by a particular behaviour of the 
spectral density: increasing the system size L for a fixed number N of channels, we end 
up to have: 

K i^i < < ■ ■ ■ < i^iv (16) 

which allowes us to simplify V2(fa, I'b), since In | cosh(i/o) — cosh(^';,)| ~ max(i/a, i/f,), to get 
approximately independent gaussian fluctuations of Ua centered at 

= (« - 1 + (17) 

and of variance 2L/{(3Nl). Since vi — 2L/^, we see that the localization length is pro- 
portional to (5: 

e = m- (18) 

The A^-dependence of ^ is characteristic of the quasi-ld limit where tranverse locahzation 
is absent. Change of ^ and change of the magnitude of the conductance fluctuations by a 
symmetry breaking effect arc the two sides of the same phenomenon, since ^ is defined from 
the ensemble average < log{g) >~ '2L/^ and not from < g >. The characteristic cross- 
over fields reducing the magnitude of the conductance fluctuations and increasing ^ in the 
absence of spin-orbit scattering are typically larger in insulators than in conductors. The 
change from orthogonal to unitary symmetry occurs when a flux quantum $o is applied 
through the phase coherent domain in a conductor (temperature dependent) and through 
the localization domain for an insulator (temperature independent). The corresponding 
fields B* are of order of fa 10 — 100 gauss in a conductor at very low temperature (a few 
mK) and 0.5 — 2 teslas for a large locahzation length insulator. For insulators with 
strong spin-orbit scattering , if one assumes that the removal of the Kramer degeneracy 
for z/i by an applied magnetic field B introduces a splitting which is negligible compared 
to the large value of 2L/^, the transition /3 = 4 — > /? = 2 essentially halves the value 
of ^. The existence of a positive magnetoresistance for insulators with strong spin-orbit 
scattering is in this theory the counterpart of the suppression of weak-anti-localization in 
conductors, in agreement with some experiments [16] and in disagreement with directed 
path theories [10]. 

The second important result that we have derived concerns the sample to sample 
fluctuations of the logarithm of the conductance since 

- \og{g) ^ 2L/i = i/i (19) 

in this quasi-ld insulating limit. We obtain a normal distribution of log(g') with 

var{log{g)) = - < log{g) > (20) 

Note that contrary to approximations based on directed paths, log{g) has a one-parameter 
distribution. 

In summary, the above results indicate that the validity of a "Coulomb gas model" 
for the {Aa} in the quasi-lo? limit is not at all equivalent to the validity of Wigner-Dyson 
statistics since the spacings between successive Aa become so exponentially large that 
the two-body interaction is reduced to an effective one body potential. In particular, 
this density effect is responsible for the non-universality of var{g) in insulators. But one 
cannot check easily in this limit the underlying presence of the logarithmic interaction. 
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If we consider now strongly disordered L* L* L cubes, the increase of their size L will 
increase the number N — {kf.L)^~^ of channels in addition to the longitudinal length. 
Real 3d-insulators are then more likely than quasi-ld systems to be characterised by an 
exponentially large value of Ai in the large L-limit without having exponentially large 
intervals between successive {\a}- Does the Coulomb gas analogy which is at least an 
excellent approximation in the absence of tranverse localization, remains of some interest 
in a strong disorder hmit? Do two successive Aa still repel each other when their spacing 
is not exponentially large, as imphed by the ansatz? Does the coulomb gas "temperature" 
take the classical values implied by symmetry considerations? This is the first issue that 
we have to check in our numerical study. 

4. Density and confining potential in the presence of transverse localization 
(3d-limit) 

We have calculated series of {Aa} for different reahzations of a microscopic model. We 
have used a 3d-network of disordered quantum wires as explained in appendix 1. Each 
link of this network is characterized by a one-dimensional uniform electrostatic 

potential which is taken at random with a rectangular distribution of width W. 
Real (propagating mode) or imaginary (decaying mode) wave-vector kij = (E — eijY^"^ 
characterizes the electronic wavefunction along the link The matching conditions at 

each intersection satisfy the basic symmetry and conservation laws of quantum mechanics. 
We have alternatively used the Anderson model (1) with the same distribution for the 
random potential e^, obtaining similar conclusions. (See also another recent numerical 
study of 3d Anderson models by Markos and Kramer [23]). 

In figure 1, the density (t(z/) is given for 3d-cubes close to the mobility edge (W=4) 
and for larger disorder (W=7). Contrary to disordered conductors where this density is 
mainly concentrated near the origin [20], one can see that for large values of W, 
is negligible in the vicinity of the origin, indicating exponentially small transmission. In 
figure 2, we show the ensemble averaged positions < Ua > which are not widely separated 
for a dielectric cube of size L, contrary to the quasi-ld limit. One can note however that 
the spacings are larger near the lower edge of the spectra than in the bulk. We note that 
g < 1 both for W = 4,1/1 ^ 2 and W = 7,Ui ^ 10 when the system size L = 6, but the 
study of the size dependence of the conductance g{L) gives j3{g) = dlog{g) / dlog{L) ^ 
for W=4.5 (mobility edge in the thermodynamic limit). 

To complete the Coulomb gas model, one needs to know the one-body potential V{\). 
In figure 3, V{\) obtained from an ensemble of numerically calculated series of {Aa} is 
given. Assuming the mean-field expression (5), we have directly performed the ensemble 
average of Yl,a=i In |A — Aa|, that we compare to an analytical fit whose the functional form 
is derived in appendix 2: 

y(A) = a.ln2(l + A). (21) 
Ao 

The parameters a and h — I/Aq are adjusted to optimal values given in figure 3. 

5. Short range repulsion of the transmission modes and sensitivity to an 
applied magnetic field. 

The Coulomb gas model implies that successive {va} interact as: 

In I cosh(i/a) — cosh(i/6)| a; In — i/a\ (22) 
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if they are large and if their spacing is smalL This logarithmic interaction must prevent 
small spacings (level repulsion) and the Wigner-surmise should give an approximate fit for 
the distribution of their spacing measured in average spacing units. We have studied this 
distribution, both for the first spacing (z/2 — ui) and for spacings between i^-levels in the 
bulk of the spectrum. As shown in figure 4, level repulsion persists even in the presence of 
transverse localization. For < ui >^ 2 the obtained histogram is very close to the usual 
prediction of random matrix theory (Wigner surmise for (3 — 1). This distribution also 
describes the spacings in the bulk of the spectrum (e. g. — uq) where the associated 
decay lengths IL/va are quite small. For larger disorder {W = 1, < vi >^ 10), level 
repulsion is weaker. The obtained histogram can then be fitted by an empirical formula 
proposed by Brody [24] (with q = 0.7) 

Pg(s) = As" exp[-Q;s«+^] (23) 

in order to interpolate from the correlated random matrix spectra (Wigner surmise q = 1) 
to uncorrelated levels (Poisson distribution q = Q). We can see in figure 2 that the 
mean spacing U2 — I'l ~ 4, which can partly explain the weakness of the repulsion 
(In I cosh(z/2) — cosh(z/i)| is more sensitive to the large value of z/2 than to the value of 
^2 ~ ^i)- However, looking deeper in the spectrum, we can still note some deviation of 
the actual distribution of (1/7 — i^q) from the Wigner surmise (figure 5), while the mean 
spacing is much smaller. But the persistence of a strong eigenvalue repulsion between 
transmission modes with extremely short characteristic decay lengths (^7 = 2L/i/7 Ri 0.5) 
is by itself quite noticeable: the spacing distribution being much closer from a Wigner 
surmise than from a Poisson distribution characterizing uncorrelated eigenvalues. 

Since we have approximately obtained spacing distributions in the presence of time 
reversal symmetry (orthogonal case) close to the Wigner surmise, it is interesting to know 
what is the value of the cross-over field B* needed for changing this distribution to the 
Wigner surmise characteristic of the unitary ensemble. It is also particularly interesting to 
see if B* characterises also the typical field scale where a change of ^ occurs, in agreement 
with our theory relating the value ^ to the importance of the fiuctuation of the {Aq} which 
are controlled by the "temperature" We have apphed increasing magnetic field B 

(measured in units of $0 per lattice cell). 

We first note that weak fields (about one fiux quantum through the sample) have no 
clear effects on the different < Ua > and on the spacing fluctuations. Near the mobility 
edge (W = 4), we observe as expected both the Wigner-surmise of the unitary ensemble 
{/3 = 2, figure 6) and a shght change of ^ (increase by a factor fa 1.3) for B — 0.25. 
One can see also in figure 6 that P{u2 — ^^i) does not depend very much on B for larger 
disorder. We think that this is partly due to the gauge invariance of a lattice model 
which makes impossible to remove the coherent backscattcring effects within too small 
localization domains. Two conclusions emerge from the magnetic field dependence of the 
spacing fiuctuations. First, the cross-over field B* is much higher in insulators than in 
metals. Contrary to a metal where the typical electronic trajectories are random walks 
exploring the whole coherent sample, the typical trajectories in an insulator are enclosing 
much smaller areas characterized by much larger dephazing fields. Our results agree with 
the criterion B* pa $o/C^- Second, the necessary field in order to significantly change 
^ is precisely the same which changes the spacing distribution, in agreement with the 
existence of a relation between the change of ^ and the change of p. 
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The study of the spacing distribution is not the more accurate test of the vahdity 
of the random matrix ansatz for 3d-insulators, probing more the short range eigenvalue 
repulsion than the long range spectral rigidity. The usual statistical measure of the long 
range rigidity is the Aa-statistics introduced by Dyson and Mehta [25]. Once the spectra 
have been numerically unfolded [26] to series of constant averaged density, this statistics 
measure how the rescaled series deviate from a regular sequence as a function of the 
average number of considered levels. We show in figure 7 how the Aa-statistics calculated 
from the {Aa} of a 3d Anderson model deviates from the random matrix behaviour when 
^ 1. Similar behaviours characterize the network of random quantum wires previously 
considered. 

The deviation of P{S) and of the As-statistics to the usual random matrix predic- 
tion can be partly due to the highly non-uniform density of the relevant variables {Aa} 
which could strongly weaken the effect of the logarithmic repulsion for the rescaled vari- 
ables. In the quasi-lrf limit for instance, we have shown that the ansatz (8) yields quasi- 
independent level fluctuations, giving a P{S) or a Aa-statistics closer to Poisson than 
to Wigner. In references 27, a crossover in P{S) and in A3 from Wigner to Poisson is 
obtained within the random matrix framework, assuming infinite range logarithmic re- 
pulsion and a one-body potential qualitatively similar to Eq. 21. However, since we 
have shown for disordered conductors [20] that the 3(i-value of the amplitude of the con- 
ductance fluctuations is contradictory with infinite range logarithmic repulsion for the 
{Aa}, wc speculate that the long range part of this repulsion must be even more severely 
damped by transverse localization than by transverse diffusion. But, from the study of 
those two microscopic models, we see that a short range repulsion persists in the presence 
of transverse localization. The shorter is C,, the shorter is this range. 

We note that a dimensionality dependent damping of the logarithmic repulsion has 
been shown [28] to be consistent with the known density-density correlations of Hamil- 
tonian spectra in disordered conductors. Similarly, the long range part of the two body 
interaction in (8) may need to be modified in order to still write the distribution P({Aa}) 
as the Gibbs distribution of a fictitious hamiltonian. 

6. Disorder dependence of the typical conductance and of its fluctuations 

In this section, we compare the actual behaviour of < ln{g) > and of its fluctuations 
to an extension of the results of section 3 and to the predictions implied by a direct path 
approximation [6-10]. We first note that contrary to the quasi-ld hmit where — In(gr) Ri ui, 
— \n{g) is slightly larger than i/i in three dimensions ( 1/2,^3,... weakly contribute to 
X]^i(cosh(z/a) — 2)~^). Extracting the localization length from the size dependence of 
< ln{g) > and locating the critical value Wc from the criterion dln[g)/dln{L) = we 
obtain (figure 8) a simple linear dependence: 

r\W) (x{W- Wc) (24) 

of the inverse localization length for the 3(i-network of quantum wires. In the vicinity of 
the critical point, this linear law is consistent with the value reported in the experimental 
litterature for the critical exponent, but lower than most of the other numerical results 

given by finite size scaling studies of various microscopic models (for a review, see Ref. 
29). But we mainly point out in this study that this linear dependence persists when the 
localization length reduces to the lattice spacing, without showing a logarithmic behaviour 
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that one could naively guess from the length dependence h of the transfer term in a 
directed path approximation (Eq. 4). 

Considering the sample to sample fluctuations of ln{g) , we have obtained histograms 
in rough agreement with a normal distribution of the log-conductance. The variance of 
ln((?) varies as a function of the disorder and of system size such that relation (20) continue 
to be observed (figure 9) in three dimensions. A similar result have been noticed in the 
2d-localized regime [21] and is consistent with one-parameter scaling, and inconsistent 
with directed path theories for the investigated range of parameters. This one-parameter 
law is observed in different microscopic models (Anderson model and network of random 
quantum wires). The comparison between the "random matrix prediction" and the di- 
rected path result can be sligthly biased since we consider the side to side transmission 
through a disorder cube, and not how electrons are transmitted from a corner to the 
opposite corner of a cube. We have noted also that the equality between the mean and 
the variance of i/i is even more accurate than for In(gi). 

7. Magnetoconductance Fluctuations and Sample to Sample Fluctuations. 

In disordered conductors, a magnetic field B.L^ ^ $o removes time-reversal symmetry 
(cross-over from the orthogonal ensemble to the unitary ensemble), suppressing weak- 
localization corrections and halving the variance of the universal conductance fluctuations. 
A (typically weak) magnetic field change of this order induces also magnetoconductance 
fluctuations of magnitude equal (to first order in 1/g) to the sample to sample fluctuations 
under applied magnetic field [4]. We have already noticed [15] in two dimensions that 
this ergodicity remains essentially correct in the localized regime for (large) field change 
AB Ri B* — $o/C^, at least when —log{g) fa 2 — 7. This "ergodic" behavior was proposed 
for explaining the large reproducible magnetoconductance fluctuations observed in the 
Mott hopping regime even for macroscopic size insulators (L ^ 1mm). It is then inter- 
esting to compare in three dimensions the typical magnitude of the magnetoconductance 
fluctuations on a given sample to the sample to sample fluctuations at given fleld. We do 
not mean here a rigorous statement about possible exact ergodicity in insulators but a 
mere comparison of the typical orders of magnitude. Due to the gauge invariance of our 
system under change of applied flux through lattice cell by $o and due to the large value 
of B* in insulators, a large statistics of those magnctofluctuations cannot be obtained. 

In Figure 10, we show how the flrst {ua} (o = 1, 2, . . . , 5) fluctuate for a given sample 
with W — 5.5 as a, function of B. The fluctuations of ui essentially agree with an ergodic 
behaviour, as illustrated in flgure 11. In addition to the large ergodic fluctuations induced 
by AB B*, one can notice smaller non ergodic fluctuations which characterizes also the 
other Pa deeper in the spectrum (figure 10) or Ui itself for larger disorder (figure 12). We 
note that larger disorder increases the scale B* of the "ergodic" fiuctuations and prevents 
us to observe them when vi ^ 1, leaving only smaller fluctuations of characteristic fleld 
scale which does not seem to depend on W. 

We point out that B* in a conductor is temperature dependent (the relevant area 
being the coherent domain), while B* in an insulator corresponds to applying a flux 
a; $0 through a temperature independent localization domain typically smaller than the 
coherence length fa Lmou in an activated regime. 

8. Effect of a variation of the Fermi energy 
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Conductance fluctuations can be observed in insulating wires where the carrier con- 
centration can be controlled by a gate potential, as shown in Ref. 18. It is interesting 
to compare the "gate" -induced fluctuations observed at very low temperature with the 
Fermi energy dependence of g occuring in our microscopic model. Moreover, this energy 
dependence provides interesting information about characteristic times of the transmis- 
sion processes. The simpler case is an isolated resonance describable by a Breit-Wigner 
formula [30]: 

where is the energy of the resonant state and the inverse of the resonance width Fe 
gives the (long) time spent by the electron trapped in this resonant localized state. In 
more complex cases where resonances and direct processes arc mixed, a smooth energy 
dependence of the transmission fluctuations still corresponds to predominant fast direct 
processes, while short energy scale fluctuations are related to slow resonant transmission, 
a priori neglected in a directed path approximation. In flgures 13. a and 13. b, one can 
see how the "gate" -induced fluctuations are different when g ^ 1 (mobility edge, roughly 
symmetric "normal" fluctuations) and when the system is deeper in the localized regime 
(a few resonant peaks superimposed on a smooth exponentially small background). Both 
in the experiments reported in Ref. 18 and here, these peaks can be fitted by the Breit- 
Wigner formula (25). Of course in the experiments, inelastic processes must in general be 
taken into account [17]. In figure 14, we can see how the applied magnetic field modifies 
this resonant behaviour. When there is a resonant peak, a variation of B detunes the 
resonance from its i? = value, inducing large conductance fluctuations around that 
peak. On the contrary, far from the resonances, in the "valleys" , the S-dependence is 
much weaker. This leads us to distinguish two effects of a variation of B. A (small) 
orbital effect which modiflcs the quantum interferences of the shortest directed paths 
which controlled a long range hop without intermediary resonances and a large effect 
resulting from the shifts of the resonant peaks. When the sample to sample fluctuations 
are studied, these two efects are mixed, though the probability to have a tuned resonance 
in a sample taken at random is very small. 

A more systematic study of the Wigner time delay matrix will be reported in a fu- 
ture publication for separating the effects of the long resonances from the directed path 
contribution. 



9. Compcirison with level statistics of the Hamiltonian matrix. 

In this study, we have investigated the validity of a (too) simple maximum entropy 
description of the transfer matrix in the limit of extreme (transverse) localization. In 
this description, the logarithmic interaction between eigenvalues and the value of j3 are 
related to a good randomization of the eigenvectors. Therefore, such a description is not 
appropriate for matrices whose eigenvector distribution undergoes important change as a 
function of the system parameters (e.g. eigenvector localization). 

For the hamiltonian matrix, the understanding which emerges from recent studies [28, 
31] can be summarized as follows in three dimensions. The logarithmic repulsion remains 
valid for a, N (x L nearest neighbour levels in the metallic phase, for a cube of size Li^. In 
the localized phase, some level repulsion is still observed for small L, but is expected to 
disappear in the thermodynamic limit where the levels would be totally uncorrelated as 
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those of hamiltonian (1) in the absence of hopping terms. One can argue that there is a 
poor overlap of the wave fuctions due to exponential localization, such that levels mainly 
depend on spatially disconnected part of the random substrate. Therefore levels are 
less and less correlated for large L in the absence of long-range correlation of the random 
variable in (1). For finite size, one has intermediary cases where the spacing distribution 
P{S) differs from the Wigner surmise more or less as in figure 4, and the deviations scale 
with L according to usual finite size scahng. The mobility edge is characterized [32] by a 
novel scale invariant P{S). 

It is then clear that, in the case of the hamiltonian matrix, "classical" random matrix 
theory with eigenvector distribution invariant under orthogonal, unitary or symplectic 
transformations and without eigenvector-eigenvalue correlations have to be extended to 
more complicated theories. Banded matrices with independent entries could be more 
realistic for describing hamiltonians in the quasi-ld limit, but very few analytical results 
are available for their eigenvalue distributions [33]. 

The random matrix X differs from the hamiltonian matrix. Localization appears in the 
presence of exponentially large eigenvalues {A^}. Our simple maximum entropy approach 
contains, via the average density cr(A), some information on the localization of the wave- 
functions, contrary to a similar approach for the hamiltonian matrix. One can then expect 
a better approximation of spectral correlation functions from this type of approach for 
X than for H. Nevertheless, a numerical study of the X-eigenvectors and of the scale 
dependence of P{S) in the localized regime will be very useful. The persistence of level 
repulsion for X is associated with the spatial overlap of electronic paths associated with 
neighboring channels. This overlap can be partly related to the importance of resonant 
tunnelling through a few localized states. Electronic modes described by the first few 
would be essentially modes where the electrons would be trapped during a long time in 
overlapping sets of certain favorable resonant states before being transmitted, as indicated 
by the large energy fiuctuations shown in section 6. Single scatterer motions in those 
sensitive regions would induce correlated fluctuations of this set of i/a, as observed in Ref. 
12. On the contrary, in a situation where shortest directed paths are predominant, one 
can imagine that the modes of transmission from on side to the opposite side of a large 
cube are associated to different non intersecting paths and are essentially uncorrelated. 

10. Compcirison with experiments and activated transport. 

A comparison between random matrix description, directed path theory, numerical 
simulations and low temperature experiments has to be done with great care. We con- 
sider in this work quantum transmission without activated process and we do not a priori 
neglect the important phenomenon of resonant tunneling which could be more or less 
truncated at finite temperature by activated hopping. We refer for this to the join paper 
by Bouchaud and Ladieu [17], where the cross-over between quantum coherent trans- 
mission and activated transport is considered. Let us mention the physical argument 
implicitcly behind the directed path approximation: if transport is not entirely coherent, 
but partly activated, one is precisely interested by long quantum hop between nearly res- 
onant states, without additional resonances between them. In other words, the lifetime of 
an elastic resonance could be too long compare to the phase-breaking time at finite tem- 
perature. Then the propagator evaluated in directed path approximations is implicitely 
related to activated transport and might not necessarily describe zero temperature elastic 
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quantum transmission. In one dimension for instance, Azbel has pointed out the role of 
resonances [34] and shown clearly that the shortest chains have not necessarily the best 
transmission ( a longer chain with a tuned resonance can be more transparent for the 
electrons than a shorter one). 

When the localization length is very short (of order of the Bohr radius of donor or ac- 
ceptor impurities), we are dealing with an "usual" bounded state and long range quantum 
transmission over Lmou is restricted to interferences between forward scattered paths in 
a "cigare shape" domain of longitudinal length Lmou ^-nd of tranverse length {Lmou-CY- 
Those directed hops are characterized by a relatively short time. When the localization 
length is longer, we are dealing with a more complicated Anderson insulator, and one 
can question if the long time elastic transmission that we study in this work can occur 
between successive inelastic collisions. 

However, despite the difficulty to take into account activated process and complicated 
quantum resonances, an element of answer is given by some striking similarities between 
the energy and magnetic field dependence of the conductance of our zero temperature 
microscopic model and the behaviours observed by Ladieu, Mailly, and Sanquer on a 
Si:GaAs wire at very low temperature, as illustrated in figure 15. 

Appendix 1: Microscopic 3d-model of disordered quantum wire networks 

The studied microscopic model can be considered as a system of narrow wave guides 
which intersect each other. We take the limit in which the transverse size of each wave 
guide tends to zero so that the motion between intersections is virtually one dimensional. 
In this limit, there is no underlying simple Hamiltonian. The mathematical foundation of 
the model has been studied by Exner and Seba [35] . The basic symmetry and conservation 
laws of quantum mechanics determine the motion completely once the matching conditions 
at each intersection arc specified. Similar concepts have been used by Avion and coworkers 
in their study of adiabatic networks [36]. For simplicity of presentation the number of 
space dimensions is fixed at d = 2 for the illustrative part (figures etc). Our numerical 
calculations are however carried out for a three dimensional network d — 3. Consider 
then a square lattice of one dimensional wires with lattice spacing a. We will denote by L 
the longitudal dimension of the system and by its transverse dimension (our numerical 
calculations are performed on cubes with L = N). On this lattice of wires an electron can 
move subject to the laws of quantum mechanics. In the most general case there is a wave 
number kj on each bond j of the lattice. We have kj = ^Ef — Vj where vj are random 
numbers drawn from a rectangular probability of width W. li W > Ef, some of the kj 
are purely imaginary and there is a disordered combination of plane-wave propagation 
and tunneling. In Fig. (16) we depict such a finite lattice. A site can therefore be labeled 
by a pair of numbers {£,n). £ = 1,2, L;n = 1,2, N. An electron can travel freely 
on each bond (wire) subject to the Schrodinger equation. We will employ the index j for 
a general bond (horizontal or perpendicular) and n for horizontal bonds only. Thus, if 
the wave number on a certain bond j (j=l,2,..LN) is kj, the electron's wave function on 
this bond is a linear combination of plane waves 

= aje'^^"" + hje-'^^'', (26) 

where for horizontal (vertical) bonds, the coordinate x is measured from the left (down) 
site. Imagine now that the A^ sites on the left of the lattice (these are the sites (1, r^)) and 
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the sites on the right of the lattice (these are the sites (L, nj) are connected to electron 
reservoirs by free wires, namely, the electron's wave number on these wires is the Fermi 
wave number kp- Then we have a quantum mechanical system whose conductance g (at 
zero temperature) can be defined in the following way. Let an incoming wave e^'^^ of unit 
amplitude reach a site (l,n) on the left column of the lattice. Then the wave function at 
the m*^ exit wire on the right (namely to the right of site (L,m)) will be tmnC**^^ where 
tmn s-re the elements of the complex transmission matrix t. Likewise, the wave function 
at the m*'* exit wire on the left (namely to the left of site (1, m)) will be r^„e~''^^ where 
rmn are the elements of the complex reflection matrix r. Unitarity (current conservation) 
implies the equality 

tt^ + rr^ = th + rV = /jy, (27) 

where J^v is the unit NxN matrix. The linear response theory now relates g to the 
transmission coefficients in a simple way. If we consider this system as a collection of 
random resistances we notice that, unlike the classical situation in which conductances 
combine rationally to one another (either in parallel or in series) the situation here is 
much more complicated, since at each intersection the wave function must be matched 
according to the rules of continuity and current conservation. Let ^i, ^2, ^3 and ^4 be 
the wave functions of the electron in the four links connected by a given site. The four 
respective wave numbers are denoted by /ci, k^i kz and k^ respectively. Then at this site, 
the following continuity equation should be satisfied, 

= ^2 = ^3 = ^4- (28) 

The relation between derivatives at the intersection is somewhat arbitrary as long as it 
is compatible with current conservation. For simplicity and consistence with previous 
works [37-38] we will assume the following relation, 

% ^% = %^% (29) 

where, following Eq. (26), 

= kjiajc'''^'' - bje-'''^^). (30) 

These relations give four equations relating the coefficients aj,bj of Eq. 26 of the plane 
waves on each of the four links (see Fig. 16). The most natural algorithm by which 
one can evaluate the conductance is the transfer matrix method. Here, however, we 
intend to study the system also in the far insulating regime. A straightforward transfer 
matrix approach is hence numerically dangerous since the product of many matrices may 
blow up. Furthermore, we intend to study a cubic system with L up to 12. For a 
3(i-cube (A^ = L) the size of the (complex) transfer matrix will then be 2N'^ = 288, 
which is quite large. It is then indispensable to invent a method which does not rely 
on a (possibly divergent) product of large matrices and at the same time will employ 
matrices of smaller size. This method is explained below in a heuristic manner. Its 
rigorous justification goes back to the theory of wave guides. It involves the algebra of 
transmission and reflection matrices whose size is N'^. Recently it has been adapted for 
studying conductances in the insulating regime [38]. The disadvantage here is that each 
step requires an inversion procedure instead of a simple multiplication as is required by 
the transfer matrix algorithm. The advantage is that the transmission and reflection 
matrices are bounded through the unitarity relation and hence there is no danger of 
spurious exponential Lyapunov divergence. 
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Consider then a wave approaching a system of two barriers "1" and "2" . The trans- 
mission and reflection amphtudes for each individual barrier i{i — 1,2), independent of 
the other one, are assumed to be given by ti and n if the wave approaches the barrier 
from the left and by t • and if the wave approaches the barrier from the right. For a N- 
channel problem, the transmission and reflection amplitudes are * A'^-complex matrices. 
They are usually arranged in an ^'-matrix pertaining to barrier i, 



Ti t' 



,i = l,2. (31) 



The 5'-matrix is unitary: Si^Si — I, S'iS'i^ = / (the unit matrix in the appropriate space). 

The transmission and reflection amplitudes pertaining to the combined system are ti2, 
t'i2 and r'12- Our task is then to express these four amplitudes in terms of the amplitudes 
of the individual channels. These relations are given explicitly as follows, 

ti2 = i2(l - r[r2)-% (32) 

ri2 = ri + t[r2{l - r[r2y^ti (33) 

t'u = t'lil - r2r[)-H'2 (34) 

r[2 = r'2 + t2r[{l-r2r[)-H'2. (35) 
In terms of the S matrices we then write 



= Si * S2, (36) 

where the star product operation deflncd through Eqs. (32 - 36) has been applied by 
Redheffer in the theory of waveguides [39] . It is then possible to add the barriers one after 
the other until the entire system is worked out. The unitarity bound on the reflection and 
transmission matrices assures the existence of the inverse (1 —r[r2)~^ and guarantees the 
convergence of the procedure for any number of steps. It is worth mentioning here that 
if the barriers are inflnitesimally close to each other, the star product is transformed into 
a non linear differential relation which eventually leads to a matrix Riccati equation for 
the reflection amplitudes. 

We have thus explained how the transmissions of individual units are combined to give 
the transmission of the whole system. Now we have to explain what are the pertinent units 
and how the transmission and reflection through each unit is computed. Transmission 
and reflection occur at each column of sites followed by the accumulation of phases as the 
wave propagates between two such columns. Thus, each unit consists of a column of sites 
together with the hnks located on its right side (see Fig. 16). In order to evaluate the 
transmission and reflection at a given column of sites we number the sites from 1 to N, 
and consider a site m(m = 1,2, ...A^). The wave number to its left is km while the wave 
number to its right is Pm- The wave numbers in the transverse direction are qm (above the 
site) and Qm-i (below the site). We adopt here periodic transverse boundary conditions 
so that go — Qn- On a transverse link between sites m and m + 1 the wave function is 
given by 

*(^/) = c„e*^'-^ + rf^e-*«-^ (37) 

where the coordinate y is measured from link m. An incoming wave of unit amplitude 
reaching site n from the left will provoke, in any site m, leftward reflected waves with 
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reflection amplitude rmn, as well as rightward transmitted waves with transmission am- 
plitude t^n- The matching equations [Eqs. (28 - 29)] now imply the following relations 
(the dependence of the transverse coefficients c and d on the incoming index n will not 
be indicated. 

where Smn is the Kronecker delta function. The equality Smn + f^mn = tmn results simply 
since reflection and transmission occur at the same point in space. From the flrst set of 
equations (39) one can now eliminate the coefficients of the transverse motions pair by 
pair using the relations 

pigma I J -igma — f fAI ^ 

We now substitute the solution of equations (40 - 41) into equation (39) and replace 
by ~ ^mn- As a result we get a set of equations for the transmission matrix t — tmn 
of the pertinent site of columns which reads 

At = 2K, (42) 

where the N * A^-matrices A and K are given by 

A{m, m)^km+Pm + i[qm cot{qma) + g^-i cot{qm-ia)] (43) 

A{m, m + 1) = -iqm csc(gma) (44) 

A{m, m-1) = -iqm-i csc(gm-ia) (45) 

K = diag{ki, k2, ...kN), (46) 

with m = 1, 2, . . . , and 1 or arc defined modulo A^. The solution of Eq. (42) yields 
the transmission matrix from left to right and then the corresponding left to left reflection 
matrix r — t — In {In is the unit matrix) is obtained. To get the matrices pertaining to 
waves approaching from the right, one has to interchange the roles of km and Pm- Since the 
matrix A is not affected by this swap, it does not require any additional matrix inversion 
procedure. It is worthwhile noticing here that in the presence of a constant magnetic fleld 
B perpendicular to the direction of propagation the only change in the above formalism 
concerns the non diagonal elements of the matrix A. If the columns of sites are numbered 
as 1, 2, . . . , L, and the magnetic flux per plaquette is (j) — Ba^ in units of $o) then at 
column ^ (^ = 1, 2, . . . , L) one has 

A{m, m + 1) ^ A{m, m + l)e^^'^'t' , (47) 

A{m, m - 1) ^ A{m, m - l)e-'^^'^'t' . (48) 

These phases are identical to those multiplying the hopping integrals (the Peierls substi- 
tution) used e.g in tight binding models. What is left in order to complete the calculations 
of the individual unit shown in Fig. (16) is to take account of the phase between two units. 
This is easily accomplished using the matrix P=diag(pi,p2, ■■■Pn) and the substitutions 

t e'^t, r^r,t'^ t'e'^, r' e'^r'e'^. (49) 
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We have thus explained how to evaluate the transmission and reflection amplitudes for 
each individual unit in the random link model. Together with the composition law ex- 
plained in the first part of this section (Eqs. (33 - 36)) we have at hand an algorithm for 
a numerical study of the full system. We repeat here that our approach will mostly be 
used in three dimensions, for which the relevant equations are slightly more complicated. 

The radial parameters {Xa} are then obtained from the eigenvalues of t.t^ using 
T„ = l/(1 + A„). 

Appendix 2. Functional form of the analytical fit for ^(A) 

Although the actual density (7(i/) differs from a simple uniform density between two 
nonzero values, it is interesting to to investigate what V{X) is yielded from this crude 
approximation (a uniform density for the z/-variable between two strictly positive values 
2L/A and 2L/B). One gets: 



/•2L/B 

Uiu) = C. / ln\ coshfz/) - cosh(u')\du' (50) 

J2L/A 



where C = {N/2L){1/B - l/A)-\ 
The derivative is given by: 

^ TT(,A ^ r 1n If exp(2V-4) -exp(;/) exp{2L/ B) - exp -{u) 
dv^ ' '^exp(2L/A) -exp-(j/)^'^ exp{2L / B) - exp{v) ^ ' 

which presents spurious divergences at the band edges 2L/A and 2L/B. Ignoring those 
divergences and after integration, one can see that the corresponding U {v) is characterized 
in presence of transverse localization (^4 7^ 0) by a flat part lor < v < vq ^ 2L/ A and 
by a quadratic increases oc (z/ — ^0 — ^- Coming back to the original variable A, 

we then expect a confining potential increasing as (ln(A) — In(Ao))^ for A S> Aq ^ 1 and 
presenting a fiat part for A < Aq. This suggests that we try to fit the V^(A) calculated from 
a microscopic 3(i-model assuming relation (5) with the following analytical expression 

y(A)=a.ln2(l + A). (52) 
Ao 

The agreement is good for small enough values of A, as shown in figure 3. It is worth to 
notice that the validity of this fit seems very general: it has been checked in 2d and "id 
disordered metals modelled either by random quantum wire networks or by tight-binding 
Anderson hamiltonians 
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Figures Captions 

Figure 1: Average density cr(z/) calculated for cubes of size L = 6 for disorder param- 
eters W = 4 (circles) and W — 7 (squares). 

Figure 2: Ensemble averaged < >, n = 1, . . . , for the same cubes as in Fig. 1. 

Figure 3: Potential V^(A) calculated for cubes of size L = 8 for different values of W 
(squares) and the analytical fit given in Eq. 21 for indicated values of a and b = Xq^. 

Figure 4: P{S) for S = {u2 — Vi)/ < z/2 — z/i > corresponding to the same cubes as in 
Fig 1-2 iox W = A (circles), W = 7 (squares) and Brody distribution with q = 0.7 (thick 
line) . 

Figure 5: P{S) for S = {ut — uq) / < — uq > for L = 6,W = 7 compared to the 
Wigner surmise. 

Figure 6 : P{S) for S — {v2 — ^i) / < 1^2 — i^i > for cubes of size L — 8 W — A (circles 
< ui >Ri 2.02) and W — 7 (Squares < ui >Ki 10.97) compared with the Wigner surmise 
for /3 = 1, 2. Apphed magnetic field B = 0.25. 

Figure 7: As-statistics as a function of the mean number n of eigenvalues A^. Unfolded 
spectra obtained in a 3(i-Anderson-model for 1^ = 30 and L = 7 at the band centre 
(^ ~ 1). The continuous curves are given by n/15 (uncorrelated spectra) and [ln(27r*n) -|- 
7 — 5/4 — 7r^/8]/7r^ (random matrix theory). 

Figure 8: The inverse localization length as a function of W — Wc for Wc = 4.5. Values 
of L beween 6 and 12 are used to evaluate ^. 

Figure 9: average (circles) and variance (diamonds) of log{g) as a function of W for 
3d cubes {L — S). 

Figure 10: The first eigenvalues oi a single sample (L = 10) at moderate disorder 
(VF = 5.5) as a function of the apphed magnetic field B. 

Figure 11: Magnetic field dependence of ui for five different samples (L = 10, 1^ = 5.5). 

Figure 12: Magnetic field dependence of ui for two samples at different values W 
(L = 10). 

Figure 13. a: Fermi energy dependence of the conductance for a single sample deep in 
the locahzed regime {g 1, L — 8,W — 6, B — 0). 

Figure 13.b: Fermi energy dependence of g near the mobihty edge {L — S,W — 6, B — 

0). 

Figure 14: Fermi energy dependence of the log-conductance of a single insulating 
sample {L — 8,W = 6) for increasing magnetic field B = 0, 0.1, 0.2, . . . , 0.5. 

Figure 15: Curve A: log{g) as a function of the gate voltage (upper coordinate) ob- 
served in a mesoscopic GaAs:Si wire (Ref.l8) where g has been multiplied by a factor 40. 
Curve B: numerical simulation of log{g) as a function of Ep using the 3(i-quantum wire 
network. 

Figure 16: Scheme of the network of quantum wires used as a microscopic model. 
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